AD/A-002  003 


A HEURISTIC  ADJACENT  EXTREME  POINT  ALGORITHM 
FOR  THE  FIXED  CHARGE  PROBLEM 

Warren  E.  Walker 

RAND  Corporation 


Prepared  for: 

Public  Health  Service 


June  1973 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


ADA00200s 


3531G2 


A HEURISTIC  ADJACENT  EXTREME  POINT 
ALGORITHM  FOR  THE  FIXED  CHARGE  PROBLEM 


Warren  E.  Walker 


June  1973 


DDCX 

i'”  DEC  IS  1974  1 1 

\:\\  ]! 

B 


P-5042 


ws'ipjmr'fcj  :;f\Tr»  ^r  A 

Approval  lor  pubHc  i***cae; 
DLs«i>tKion  Unlimitod 


Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

US  Depertmer.t  of  Commerce 
Springfield,  VA.  221  Si 


fill's  Sectfoa 

r ? E';,‘  LI 


‘■'lJ3’!0.'!,/*M!U3IUTy  COOS 
. znjir 


Any  views  expressed  in  this  paper  are  those  of  the  author.  They 
should  not  be  interpreted  as  reflecting  the  views  of  The  New  York 
City-Rand  Institute  or  the  official  opinion  or  policy  of  the  City  of 
New  York.  Papers  are  reproduced  by  The  Rand  Corporation  as  a 
courtesy  to  members  of  its  staff. 


A HEURISTIC  ADJACENT  EXTREME  POINT 
ALGORITHM  FOR  THE  FIXED  CHARGE  PROBLEM 


Warren  E.  Walker 


The  New  York  City-Rand  Institute 
545  Madison  Avenue 
New  York,  N.  Y.  10022 


ABSTRACT 


An  algorithm  with  three  variations  is  presented  for  the  approximate 
solution  of  fixed  charge  problems.  Computational  experience  shows  it 
to  be  extremely  fast  and  to  yield  very  good  solutions. 

The  basic  approach  in  all  three  variants  of  the  algorithm  is  (1) 
to  obtain  a local  optimum  by  using  the  simplex  method  with  a modifica- 
tion of  the  rule  for  selection  of  the  variable  to  enter  the  basic  so- 
lution, and  (2)  once  at  a local  optimum  to  search  for  a better  extreme 
point  by  jumping  over  adjacent  extreme  points  to  resume  iterating  two 
or  three  extreme  points  away. 

Problems  in  which  economies  of  scale  give  rise  to  separable  piece- 
wise-linear  concave  objective  functions  are  shown  to  be  easily  formu- 
lated as  fixed  charge  problems. 

The  algorithm  is  currently  being  used  by  the  U.  S.  Environmental 
Protection  Agency's  Office  of  Solid  Waste  Management  Programs  to  solve 
a problem  of  regional  solid  waste  planning:  the  selection  of  disposal  sites 

to  be  developed  and  the  determination  of  how  the  wastes  of  each  municipality 
in  a region  should  be  distributed  among  the  sites. 


Any  views  expressed  in  t'lis  paper  are  those  of  the  author.  They 
should  not  be  interpreted  as  Teflecting  the  views  of  The  New  York 
City-Rand  Institute  or  the  official  opinion  or  policy  of  the  City  of 
New  York.  Papers  are  reproduced  by  The  Rand  Corporation  as  a courtesy 
to  members  of  its  staff. 

This  work  was  supported,  in  part,  by  a grant  from  the  U.S.  Public 
Health  Service  (ES-00098). 
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IfTiRODUCTION 


The  fixed  charge  problem  can  be  stated  as 


min  z = l f (x  ) 
j-1  3 3 


Ax  = b 
x > 0 


(1) 


where  f (x^ ) 


c.x.  + k 6 

, J J j J 

JO  if  Xj=0 

(1  if  x >0 


This  type  of  problem  arises  in  many  practical  settings.  Two  of  the 
most  common  of  these  are  in  warehouse  or  plant  location,  where  there  is 
a charge  associated  with  opening  the  facility,  and  transportation  prob- 
lems, where  there  are  fixed  charges  for  transporting  any  goods  between 
supply  points  and  demand  points. 

If  all  the  fixed  charges  were  zero,  then  problem  (1)  would  be  a 
linear  programming  problem.  If  some  or  all  of  the  fixed  charges  are 
positive,  the  objective  function  z is  concave  [11]  and  it  is  well  known 
that  a concave  function,  defined  over  a convex  polyhedron, takes  on  its 
minimum  at  an  extreme  point.  The  methods  we  present  for  the  solution 
of  (1),  therefore,  will  examine  only  extreme  point  solutions. 

The  fixed  charge  problem  can  be  written  as  a mixed-integer  linear 
program  [10,  p.  253].  However,  most  mixed  integer  programming  algorithms 
are  not  computationally  efficient  for  large  problems.  This  is  true  also 
for  exact  algorithms  developed  specifically  to  solve  the  fixed  charge 
problem.  For  example,  Steinberg's  branch  and  bound  algorithm  [15] 
requires  as  much  as  47  minutes  on  an  IBM  360/50  to  solve  a 15x30  prob- 
lem. Gray’s  decomposition  approach  [8]  requires  an  average  of  lb 
minutes  to  solve  a 5x7  fixed  charge  transportation  problem  and  as  much 
as  22  minutes  to  solve  a 30-site  warehouse  location  problem  on  the  Bur- 
roughs B-5500  [9].  Murty  [13]  has;  developed  an  exact  algorithm  which 
solve!  ! with  all  k.=0  to  bound  the  total  ec>  anti  then  searches  ay;- 
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tematically  among  the  extreme  points  adjacent  to  the  LP  optimum  for  the 
minimum  total  cost.  Murty  presents  only  one  sample  problem,  solved  by 
hand.  But  Gray  [8]  tried  to  solve  one  6x8  and  six  5x7  test  problems. 

He  was  able  to  solve  only  two  of  the  problems,  running  cut  of  computer 
storage  capacity  before  solving  each  of  the  other  five. 

Since  the  currently  available  exact  algorithms  generally  require 
long  computation  times  and  large  amounts  of  storage,  a good  deal  of  ef- 
fort has  been  devoted  to  finding  approximate  solutions  to  fixed  charge 
problems.  The  fixed  charge  transportation  problem  (where  the  A-matrix 
is  in  the  form  of  a transportation  matrix)  has  been  investigated  by 
Kuhn  and  Baumol  [12]  and  Balinski  [2].  Kuhn  and  Baumol  suggest  that 
an  approximate  solution  t >,  the  problem  may  be  obtained  by  forcing  a 
highly  degenerate  solution.  This  is  accomplished  by  making  small  ad- 
justments to  the  right  hand  side  (demands  and  supplies).  The  approxi- 
mation is  a rough  one,  but  the  computation  is  quite  simple. 

Balinski  replaces  thp  non-linear  fixed  charge  objective  function 
by  an  aoproximate  linear  objective  function,  and  solves  the  resulting 
problem  using  the  standard  transportation  algorithm.  He  also  finds 
bounds  on  the  optimal  exact  solution.  This  approach  yields  a rather 
rough  approximation,  and  does  not  work  well  in  many  cases. 

Cooper  and  Drebes  [5],  Steinberg  [15],  and  Denzler  [7]  have  developed 
approximate  heuristic  adjacent  extreme  point  algorithms  for  the  general 
fixed  charge  problem.  Steinberg  and  Denzler  modify  the  linear  program- 
ming criterion  for  a vector  to  enter  the  basis,  a technique  also  used 
in  this  paper.  Their  algorithms  will  be  discussed  more  fully  below. 

Cooper  and  Drebes  modify  the  objective  function  at  certain  stages  in 
their  algorithms,  and  also  change  the  criteria  for  vectors  to  enter  and 
leave  the  basis.  At  certain  times  in  their  calculations  a vector  is  cho- 
sen to  enter  the  basis  with  the  least  fixed  charge  of  the  non-basic 
valuables.  At  other  times  a vector  is  chosen  to  leave  with  the  highest 
fixed  charge  in  the  basic  set. 

The  computational  experience  reported  for  these  methods  indicates 
that  they  will  yield  the  optimal  solution  a high  percentage  of  the  time 
and,  when  not  optimal,  they  provide  a good  approximation. 

This  paper  will  describe  another  adjacent  extreme  point  algorithm 
(with  variations)  which  appears  to  be  faster  than  those  of  [5],  [7]  and 


[15],  and  which  yields  the  optimal  solution  a higher  percentage  of  the  time. 
The  algorithm  currently  is  being  used  on  a production  basis  to  solve  large 
fixed  charge  problems  [1J,  [4].  The  algorithm  will  be  referred  to  by  the 
name  SWIFT,  for  Simplex  With  Forcing  Trials. 

THE  SWIFT  ALGORITHM 

Before  proceeding  with  cne  algorithm,  we  make  the  following  notational 
definitions: 

Xg  5 vector  of  basis  variables 
Cg  = vector  of  prices  for  xg 
kg  = vector  of  fixed  charges  for  xg 
aj  = jth  column  of  A 
B = basis  matrix 

y^  = B ^a^  (the  representation  of  a^  in  terms  of  B) 

zj H Vr 

The  algorithm  consists  of  two  sequential  phases. 


Phase  1 


The  first  phase  is  identical  to  the  standard  simplex  procedure  except 

that  the  rule  for  selecting  the  column  to  enter  the  basis  is  modified. 

In  a linear  program,  the  objective  function  will  decrease  if  the  entering 

vector,  Xj , is  selected  so  that  zj  ~ cj  > Because  of  the  fixed  charges, 

this  criterion  will  not  insure  a local  improvement  in  the  fixed  charge 

objective  function.  However,  the  non-negativity  of  a similar  quantity 

involving  - c^  and  the  fixed  charges  will  produce  an  improvement. 

Suppose  that  x^  is  to  enter  the  basis  on  the  next  iteration.  Then, 

the  leaving  vector  is  xD  , where  x„  is  determined,  as  in  the  ordinary 

Br  Br 

simplex  algorithm,  by 


min 

k 


lykj 


ykj  > 0 


*B 


(2) 


and  6j  is  the  value  which  x^  assumes  upon  entering  the  basis. 
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If  0j  > 0,  as  a result  of  such  a basis  change  the  objective  function 

i3  increased  by  k. , decreased  by  kD  , and  increased  or  decreased  by 
J Br 

9(Zj  - Cj)  depending  upon  the  sign  of  z^  - c^.  (If  0^  = 0,  the  objective 

function  remains  the  same.)  In  addition,  if  the  choice  of  xn  was  not 

Br 

unique,  one  or  more  of  the  basic  variables  which  were  positive  will 
become  zero.  In  this  case,  even  though  they  remain  in  the  basis,  the 
objective  function  is  reduced  by  their  fixed  costs.  Conversely,  it  is 
possible  that,  in  the  course  of  bringing  x^  into  the  basis,  some  basic 
variables  which  were  at  a zero  level  will  become  positive.  If  this 
occurs,  their  fixed  costs  must  be  added  in  to  determine  the  new  value 
of  the  objective  function.  This  requirement  was  neglected  by  Denzler 
and  glossed  over  by  Steinberg. 


Let 


S c 


\ \ 
yU  V 


T = 


1'\ ' °*  >u  K 0 


Then,  the  entering  vector,  x ^ , should  be  chosen  such  that 


“J  -B,  “j<Zj  “ Cj)  ' ki+  L “i 


- k^  - 0, 
r 


k,  < 0. 


(3) 


ieT 


It  is  possible  to  continue  iterating  using  criterion  (3)  to  choose  a 
vector  to  enter  the  basis  until  A^  _>  0 for  all  non-basic  columns  j. 
However,  because  the  objective  function  is  concave,  it  is  not  true 
that  when  all  A^  _>  0 a global  minimum  has  been  reached.  Even  though 
no  adjacent  extreme  point  will  yield  a smaller  value  of  z,  it  is  still 
possible  that  some  other  extreme  point  of  the  convex  set  will  be  better 
(see  [12],  p.  13). 

This  difficulty  leads  to  phase  2 of  the  algorithm — a search  for  a 
better  extreme  point  non-adjacent  to  the  current  point.  Three  closely 
related  methods  have  been  developed  for  phase  2.  Taken  together  with 


fltejpllMiiMiMrl  f 


-5- 


I 

* 


phase  1,  which  i.i  the  serve,  for  each  of  these  methods,  they  constitute 
three  heuristic  algorithms  for  the  fixed  charge  problem,  which  will  be 
called  SWIFT-1,  SWIFT-2  and  SWIFT-3.  Aside  from  the  treatment  of  degen- 
eracy in  phase  1 noted  above,  the  phase  2 procedures  distinguish  the 
SWIFT  algorithms  from  Steinberg's  and  Denzler's.  The  SWIFT  algorithms 
constitute  a more  deliberate  search  for  improvement  and  have  produced 
better  results. 

Phase  2 

At  the  end  of  phase  1,  A^  _>  0 for  all  non-basic  columns  j.  The 
phase  1 solution  may  or  may  not  be  the  optimum.  In  phase  2 one  or  more 
vectors  will  be  forced  into  the  basis,  increasing  the  objective  function, 
because  of  the  possibility  that,  by  continuing  iterations  from  a new 
polit,  the  algorithm  might  move  away  from  the  old  local  optimum  to  an 
improved  new  one.  That  is,  from  a local  optimum  an  investigation  is 
made  of  nearby  extreme  points  with  larger  objective  values  which  might 
be  adjacent  to  points  with  smaller  objective  values. 

The  three  different  methods  for  phase  2 presented  below  differ  in 
(1)  the  number  of  vectors  forced  into  the  basis  at  a time  and  (2)  the 
action  taken  if  a forcing  attempt  fails  to  produce  a better  solution. 


THREE  SWIFT  ALGORITHMS 

SWIFT- 1 (single  forcing,  non-return) 

0.  Find  an  initial  feasible  solution  to  (1). 

1.  Iterate  with  the  simplex  method,  using  criterion  (3)  to  choose 
a vector  to  enter  the  basis,  until  A^  _>  0 for  all  non-basic 
columns  j . 

a.  Let  Xq  be  this  phase  1 solution. 

b.  Let  Zq  be  the  corresponding  value  of  the  objective  function. 

2.  Force  a currently  non-basic  variable,  not  yet  tried,  into  the 
basis,  yielding  a new  solution  x^  with  objective  value  z^  z^. 

If  all  non-basic  variables  in  solution  Xq  have  been  tried  without 
an  improvement,  STOP  and  call  x^  the  (approximate)  solution, 
otherwise  go  to  step  3. 
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Iterate  as  in  step  1 until  A^  _>  0 for  all  non-basic  columns  j 
and  a local  optimum  xg  is  found. 

a.  If  x = x-  (i.e.,  no  iterating  was  possible),  return  to 
solution  Xq.  Go  to  step  2. 

b.  If  z_  < zrt  a better  solution  has  been  found.  Rename  this 
Go  to  step  2. 


s *0 
solution  x, 


c. 


0 

If  zg  _>  Zq,  go  to  step  2. 


SWIFT- 2 (single  forcing,  return) 

Same  as  SWIFT-1  except  change  step  3(c)  to  read: 

c.  If  z _>  zn  return  to  solution  x.  (the  best  solution  so  far). 
3 U U 

Go  to  step  2. 


SWIFT- 3 (double  forcing,  return) 

Same  as  SWIFT-2  except  change  step  2 to  read: 

2.  Force  an  untried  pair  of  non-basic  variables  from  solution  x^ 
into  the  basis,  yielding  a new  solution  x^  with  objective  value 
Zj  Zg.  If  all  pairs  of  non-basic  variables  in  solution  x^ 
have  been  tried  without  an  improvement,  STOP  and  call  x^  the 
solution. 


STEEPEST  DESCENT 


The  criterion  used  to  determine  the  vector  to  enter  the  basis  in 
phase  1 could  be  changed  to  a steepest  descent  criterion.  That  is, 
choose  vector  x^  to  enter  if  A^  <0  and 


A,  ■ min  A. 
3 il^B  1 


where  A^  is  given  by  (3). 

The  steepest  descent  criterion  is  rarely  used  in  solving  a normal 
linear  programming  problem  because  it  involves  finding  6^  for  each  non- 
basic  column  x^  for  which  z^  - c^  > 0.  However,  the  algorithm  described 
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above  for  the  fixed  charge  problem  requires  the  calculation  of  0^  for  at 
least  a subset  of  the  non-basic  columns.  As  a result,  this  criterion 
is  easy  to  implement  'and  adds  little  time  to  the  calculations  per  itera- 
tion, but  reduces  the  number  of  iterations  per  problem  by  between  15  and 
30  percent. 

TESi  RESULTS  AND  COMPARISONS  WITH  OTHER  ALGORITHMS 

In  order  to  test  their  heuristic  fixed  charge  algorithms.  Cooper 
and  Drebes  [5]  randomly  generated  290- (5x10)  problems  with  the  following 
properties: 


1 < c1  < 20 
1 < k < 999 

The  average  density  of  A is  50  percent. 

The  optimal  solutions  to  these  problems  were  obtained  by  complete  enumera- 
tion. These  problems  and  their  solutions  are  inc]uded  in  Steinberg's 
thesis  [16].  They  have  beer  used  by  Cooper  and  Drebes,  Steinberg  and 
Denzler  to  test  their  respective  algorithms. 

Cooper  and  Drebes  applied  their  algorithms  to  a set  of  253  of  these 
problems.  Of  these,  240  were  solved  optimally  by  their  algorithm  MI,  and 
245  by  their  algorithm  Mil.  Der.zler  obtained  optimal  solutions  to  169  out 
of  200  problems  using  his  M-l  algorithm  and  all  200  using  his  M-3  algorithm. 
Steinberg  obtained  optimal  solutions  to  235  out  of  250  problems  using  his 
Heuristic  One  algorithm,  and  255  out  of  268  problems  using  his  Heuristic  Two. 
The  SWIFT  algorithms  were  tested  on  the  22  problems  for  which  Steinberg  got 
suboptimal  solutions  using  his  algorithms.  All  three  of  the  algorithms  ob- 
tained optimal  solutions  to  the  22  problems.  In  16  of  the  22  cases,  opti- 
mality was  attained  by  the  end  of  phase  1. 
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Subsequently,  the  algorithms  were  tried  on  a random  sample  of  30 
problems  solved  successfully  by  Steinberg.  Again,  the  optimal  solution 
was  obtained  for  all  the  problems.  The  optimum  was  reached  by  the  end 
of  phase  1 in  26  out  of  the  30  problems. 

The  problems  were  solved  with  and  without  steepest  descent.  As  an 
example  of  the  effect  of  steepest  descent,  method  2 averaged  a total  of 
18  simplex  iterations /problem  without  and  15  iterations/problem  with 
steepest  descent.  Steepest  descent  also  reduced  the  average  number  of 
iterations  needed  to  reach  the  optimal  solution  by  about  3 iterations 
(in  method  2,  from  10.5  iterations  to  7 iterations). 

Table  1 gives  the  results  of  the  SWIFT  algorithms  compared  with 
those  of  Cooper  and  Drebe9,  Denzler  and  Steinberg  for  the  5x10  test 
problems. 

Cooper  and  Drebes  constructed  15x30  test  problems  by  aggregating 
sets  of  three  5x10  problems.  In  these  the  A-matrix  was  formed  as  fol- 
lows : 


where  A^,  Aj,  A^  are  5x10  matrices.  The  optimal  objective  value  for  any 
of  the  15x30  problems  is  the  sum  of  the  optimal  objective  values  from  the 
5x10  problems  associated  with  A^,  and  Ay 

Six  such  15x30  problems  were  constructed  for  SWIFT  testing.  SWIFT 
obtained  optimal  solutions  for  all  six  problems.  Optimal  phase  1 solu- 
tions were  obtained  in  two  cases.  These  two  had  A-matrices  which  were 
comprised  of  submatrices  which  had  led  to  optimal  phase  1 solutions 
to  the  constituent  5x10  problems.  The  other  four  A-matrices  contained 
at  least  one  submatrix  from  a 5x10  problem  which  did  not  produce  an 
optimal  phase  1 solution. 


..  -- nr  i miiippwpBipi  | , ipm  rr 


Table  1 

COMPARISON  OF  ALGORITHMS 
5x10  PROBLEMS 


No.  No.  X per 
Algorithm  tried  optimal  opt.  problem 


Average  no. 

iterations  Average  time 
per  per  problem 

problem  (sec.) 


Computer 


SWIFT-1 

52 

52 

100 

23 

5 

CDC  1604 

SWIFT- l-s.d.a 

49 

49 

100 

17 

0.5 

IBM  360/65 

SWIFT- 2 

52 

52 

100 

18 

4 

CDC  1604 

SWIFT-2-s .d.a 

49 

49 

100 

15 

.47 

IBM  360/65 

SWIFT- 3 

52 

52 

100 

36 

6 

CDC  1604 

MI 

253 

240 

95 

75 

20 

IBM  7072 

Mil 

253 

245 

97 

100 

20 

IBM  7072 

Both 

253 

250 

99 

175 

20 

IBM  7072 

Heuristic  1 

250 

235 

94 

6.8 

.7 

IBM  360/50 

Heuristic  2 

268 

255 

95 

15.2 

1.5 

IBM  360/50 

M-l 

200 

169 

85 

n.  a. 

1.6 

IBM  7072 

M-2 

200 

196 

98 

n • a • 

7.0 

IBM  7072 

M-3 

200 

200 

100 

n.a. 

14.0 

IBM  7072 

3 

s.d.  = steepest  descent. 

^The  add  times  (in  micro-seconds)  for  each  of  the  computers  used 
IBM  7072  12.0 

CDC  160A  7 . 2 

IBM  360/50  4.0 

IBM  360/65  1.3 

n.a.  = not  available. 
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The  15x30  problems  were  substantially  harder  to  solve  than  the  5x10 
problems.  SWIFT-2  went  from  an  average  of  15  iterations/problem  to  an 
average  of  86  iterations /problem.  SWIFT-3  required  as  many  as  392  itera- 
tions to  solve  one  15x30  problem.  However,  it  had  reached  the  optimal 
solution  by  iteration  35.  The  other  iterations  were  spent  searching  for 
a better  solution.  A comparison  of  the  results  of  using  different  algo- 
rithms for  solving  the  15x30  problems  is  given  in  Table  2. 


Table  2 

COMPARISON  OF  ALGORITHMS 
15x30  PROBLEMS 


Author 

Algorithm 

No. 

tried 

No.  % 

optimal  opt. 

Average  no. 

iterations 

per 

problem 

Average  time 
per  problem 
(sec.) 

Computer 

Walker 

SWIFT-2 

5 

5 

100 

86 

18 

CDC  1604 

SWIFT-3 

2 

2 

100 

370 

60 

CDC  1604 

SWIFT-3-s.d.a 

1 

1 

100 

376 

25 

IBM  360/65 

Cooper  & 
Drebes 

Both  MI  £. 
Mil 

70 

63 

90 

1200 

900 

IBM  7072 

Steinberg 

Heuristic  1 

90 

75 

83 

12 

2.16 

IBM  360/50 

Heuristic  2 

84 

74 

88 

43 

7.74 

IBM  360/50 

Denzler 

M-l 

22 

14 

64 

n.a. 

n.a. 

M-3 

22 

17 

77 

n.a. 

75 

IBM  7072 

s.d.  = steepest  descent, 
n.a.  = not  available. 


IPMM'*™- 
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Gray  [9]  lists  12  fixed-charge  transportation  problems  which  he  solved 
using  his  exact  algorithm.  SWIFT-2  was  applied  to  these  same  problems. 
Table  3 contains  a comparison  of  Gray's  computation  times  with  those 
achieved  by  SWIFT-2.  The  IBM  360/65  has  an  add-time  which  is  three  times 
faster  than  that  of  the  Burroughs  B-5500.  Even  correcting  for  this,  the 
SWIFT  algorithm  is  from  2 to  835  times  faster  than  the  Gray  algorithm. 

It  failed  to  get  the  optimal  solution  to  two  problems,  but  in  these  two 
cases  its  solution  was  greater  than  the  minimum  cost  solution  by  only 
.8  percent  and  1.8  percent. 


Table  3 

FIXED-CHARGE  TRANSPORTATION  PROBLEM  SOLUTIONS: 
A COMPARISON  OF  THE  GRAY  AND  WALKER  ALGORITHMS 


Problem 

no. 

Sice  of 
A-matrix 

Gray 

solution 

Time  (sec.) 
B-5500 

Walker 

solution 

(SWIFT-2) 

Time  (sec.) 
IBM  360/65 

1 

7x12 

329 

7.7 

329 

1.12 

la* 

7x12 

429 

7.6 

429 

.42 

lbb 

7x12 

579 

7.8 

579 

.43 

2 

10x24 

202 

32.6 

202 

1.36 

3 

10x24 

1999 

26.3 

1999 

1.70 

4 

12x32 

273 

171.4 

-73 

3.69 

5 

12x35 

245 

263.8 

247 

3.68 

6 

12x35 

317 

146.9 

317 

6.11 

7 

12x35 

1638 

97.0 

1668 

3.01 

8 

12x35 

2289 

3262.8 

2289 

3.89 

9 

14x48 

314 

1510.1 

314 

7.61 

9ac 

14x48 

2357 

71.4 

2357 

5.89 

identical  to  problem  1 with  20  added  to  each  fixed  charge. 
^Identical  to  problem  1 with  50  added  to  each  fixed  charge, 
identical  to  problem  9 with  250  added  to  each  fixed  charge. 


^ vtjerxrvzR*  --**.-*■  v-y.  't'^:-- 
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THE  SWIFT  ALGORITHM  APPLIED  TO  REGIONAL 


SOLID  WASTE  DISPOSAL  PLANNING 


In  1968  Skelly  [14]  developed  a model  for  use  In  planning  for 
regional  solid  waste  disposal.  The  problem,  viewed  as  a network  model, 
was  to  determine  the  disposal  facilities  (sinks  and  trans-shipment 
points)  which  should  be  developed  to  handle  the  refuse  generated  by 
several  communities  (sources)  such  that  the  total  costs  of  transporta- 
tion, treatment  and  disposal  are  minimized.  The  formulation  is 
basically  a capacitated  warehouse  location  problem.  The  cost  func- 
tion includes  linear  transportation  costs,  fixed  costs  for  using  a 
site  and  piece-wise  linear  concave  site  operating  costs  (which  can 
be  replaced  by  fixed-charge  equivalents  as  is  shown  in  the  Appendix) . 

The  SWIFT-2  algorithm  with  steepest  descent  was  used  to  solve 
the  problem  using  fifteen  different  sets  of  actual  data.  The  largest 
problem,  a 25-city,  9-site  problem  with  42  constraints  and  282  columns 
(including  slacks  and  artificials)  was  solved  in  8.9  minutes  on  an 
IBM  360/65.  Of  the  fifteen  problems,  five  found  improved  solutions 
in  the  forcing  phase,  the  rest  reached  their  final  solutions  in 
phase  1. 

This  model  has  been  further  developed  by  Roy  F.  Westlri,  Inc. 
as  part  of  a comprehensive  state-wide  solid  waste  management  study 
for  the  New  York  State  Department  of  Environmental  Conservation  [1] 
and,  more  recently,  by  the  Federal  Environmental  Protection  Agency's 
Office  of  Solid  Waste  Management  P-ograms  (OSWMP)  for  use  by  regional 
planning  authorities  throughout  the  United  States.  The  computer 
program  being  used  has  a matrix  generator  to  simplify  data  input, 
and  a report  generator  for  presenting  the  results  in  a clear  and 
meaningful  manner.  But  it  uses  the  SWIFT  algorithm  to  solve  the 
fixed-charge  problem.  OSWMP  is  currently  using  the  model  to  develop 
regional  disposal  plans  for  the  Seattle  area  and  for  an  11-county 
area  of  Texas  which  includes  the  city  of  Dallas.  The  Texas  problem 
has  200  refuse  sources,  70  potential  disposal  sites  and  800  source- 
site  transportation  pairs.  This  results  in  a constraint  matrix 
having  400  rows  and  870  columns  (not  including  slack  and  artificial 
variables) . 
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Appendix 

TRANSFORMING  A CONCAVE-COST  LINEAR  PROGRAMMING 
PROBLEM  INTO  A FIXED-CHARGE  PROBLEM 


Economies  of  scale  lead  to  cost  curves  whose  slope  decreases  as 
the  independent  variable  increases.  Because  of  the  frequency  of  its 
occurrence  in  applications , the  minimization  of  concave  objective 
functions  subject  to  linear  constraints  has  received  considerable 
attention  in  the  literature  [6,  p.  543],  [10,  Chapter  3],  ana 
[3,  Chapter  Xj.  The  proposed  algorithms  for  this  class  of  problems 
usually  involve  integer  programming  or  lengthy  search  procedures. 

However,  it  will  be  shown  that,  if  the  function  is  piece-wise 
linear  and  separable,  or  can  be  approximated  by  a piece-wise  linear 
separable  concave  functional,  the  problem  can  be  formulated  as  a 
fixed-charge  problem.  For  this  development,  we  first  assume  that 
any  strictly  concave  objective  function  has  already  been  approximated 
by  a piece-wise  linear  separable  concave  functional.  Then,  the  con- 
cave-cost linear  programming  problem  can  be  stated  as: 


min  z - l <f>  (x  ) 
j-1  3 3 


s.t. 


n 

I 

1-1 


aij  Xj  " bi 


i - 1.  2, 


(4) 


Xj  21  0 j * 1,  2,  ...»  n 


where  each  of  the  ^ (x^ ) is  piece-wise  linear 
and  concave. 


We  will  construct  a fixed-charge  problem  which  is  equivalent  to  (4). 

The  resulting  problem  can  then  be  solved  by  any  fixed-charge  algorithm. 


Suppose  that  <f>j  (Xj)  is  composed  of  linear  segments.  Let 

= slope  of  ith  segment  of  (j)^  (x^): 

CU  > C2j  > •••  > 

f^j  = y-intercept  of  ith  segment  of  <J>j  (x^) 
when  extended  to  the  y-axls: 

0 - flj  < f2j  < ‘ " < f rj j 

h^j,  j = break-points  of  <J)j  (Xj)  on  the  x^  axis; 

J h0j  5 °- 


Using  this  notation  l}^  (x^)  can  be  represented  graphically  by  the  solid 
curve  in  Fig.  1 below. 


~ ^^frijaww'S^F.  *!^m,WT!r?< ••  r***y*' "*?• J-v  *y,  v;;-- ♦*--,»  •*  .*^r,w •*  ’ 
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For  each  activity  Xj,  (j  - 1,  2 n) , define  fj  new  variables 

A , A , ....  A . Associate  with  A±j  the  variable  cost  c^  and  the 
fixed  cLt  f . ^Each  variable  A^,  therefore,  has  a fixed  cost  objec- 
tive function  of  the  form: 


Figure  1 can  be  viewed  as  the  superposition  of  these  r^  fixed  cost  objec- 
tive functions  onto  a single  set  of  axes.  The  function  ^ (x^)  is  the 
lower  envelope  of  these  cost  functions.  That  is,  for  a given  value  of  x^ , 
the  value  of  the  function  <}) j (x^)  is 

min  (c..x  + 
l<i<T.  J J J J 


where  6 


j 


j 0 if  Xj  - 0 
j 1 if  Xj  > 0 


| 


Let  Xj  - Afj  + ...  + Ar^. 
Then  (A)  can  be  rewritten  as: 


ii 

min  z,  - I “in  (c  A1;)  + 
j-1  1<  i <r  j 

s.t.  V a.  . T A..  ■ b,_  (k  ■ 1,  2 , . . . , ■) 


I l Aij  " bk 

j-1  i-1  J 


(5) 


hi-l,j5ij  5 " hij6iJ 


\ 0 if  Aaj  - 0 


ii  1 if  A,,  > o. 
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+ could 


'ij“ij  ' ij  ij 


n rj 

The  fixed-charge  objective  function  £ I (c-t  *Aj 

j-1  i“l 

be  substituted  for  the  objective  function  in  (5)  if  at  most  one  A^  would 
be  positive  for  each,  j in  any  optimal  solution. 

Consider  problem  (5)  without  the  constant  3-  Aij  hij(Sij  81111 

with  the  objective  changed  to 


n rj 

min  z - l l (c  .A  . + f 
j“l  i=*l  J J 


ij'V* 


The  resulting  problem  is: 


(6) 


We  state  and  prove  two  theorems  which  together  show  that  problem  (6) , 
which  is  a fixed-charge  problem,  is  equivalent  to  problem  (4) . 


THEOREM  1 

In  an  optimal  solution  to  (6),  at  most  one  A^  will  be  positive 
for  each  j . 

proof;  Suppose  an  optimal  solution  has  more  than  one  A^  positive  for 
some  j.  Consider  any  two  of  them,  say  A^  > 0 and  A^  > 0.  Without 
loss  of  generality,  let  a <b.  The  cost  of  this  solution  is: 


* K + (taj  + + £bJ  + °blV 


l . I < Vlp  + W + ? ■ ■ <cii*13  + Vu’ 


p^j  i"l 


i^a.b 


ij  ij  ij  ij'* 


where  K 
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Con  side  r reducing  A to  0 and  Increasing  A,  by  A . , while  leaving  all 
Hj  Dj  fij 

other  variables  unchanged.  This  Is  also  a solution,  and  Its  cost,  In 
terms  of  the  original  variables.  Is: 


Then, 


z2  ’ K + £bj  + (4«j  + V' 

z2  " Z1  ‘ 4aj  <cbJ  " CaJ*  " £aj. 


But,  by  concavity,  c^  > c^^ , and,  by  assumption,  f^  > 0 and  A^  > 0. 
Thus,  < z^,  contradicting  the  assumption  that  we  had  an  optimal 
solution. 


THEOREM  2 


If  a A^  la  positive  In  an  optimal  solution  to  (6) , Its  value 

will  always  fall  between  h.  . and  h . 

lj 

Proof:  Suppose  that  A > 0 In  some  solution  to  (6) . By  Theorem  1, 

aJ 

x.  - A ..  Ifh  , < A . < h the  theorem  Is  proved,  so  assume 

j a-l,j  - aj  - a,j’ 

that  (a)  ^ j or  (b) 

Case  (a) : A^  > h^ : The  objective  function  corresponding  to  this 

solution  Is 


- \n  \ + faj  + c.j  ‘j 


(7) 


where 


l Ic. 

1-1 


A,.  + 


Suppose  A . , . were  Increased  from  0 to  x.  and  A . were  decreased  from 

«+l»j  j »j 

xj  “ °- 

The  objective  function  for  this  new  solution  would  be 


s£££2s£2ii-:rr-; * ~ ' ™**»*m>.  iij#ui»,ij»jR,#i'i>Mi‘i.  t juu,h 
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z .,  ■ V k + f .,  . + c ,,  . x. . 
»fl  a afl.j  «fl,j  j 


(8) 


Subtracting  (7)  from  (8)  produces 


S*fl  “ Za  ” (fa+l,j  " faj'  + *Ca+l,j  " Caj}  * j * (9) 


At  break-point  h^,  the  contributions  to  the  objective  value  from  A^  and 
A^+1  are  the  same;  that  is 


fij  + Cij  hij  * fi+l,J  + Ci+l,j  hij 


or 


(fi+l,j  ‘ fi,j)  “ (cij  ” Ci+l,j)  hij * 


(10) 


Substituting  (10)  into  (9)  we  obtain 


za+l  ” za  “ (caj  “ Ca+l,j)  (haj  " V* 


By  concavity,  caj  > c^1#J  - - c^1#j 


> 0. 


By  assumption,  h , < x.  h . - x.  < 0. 

* j j ® » J J 


Thus,  z < z , contradict.’-  g the  assumption  of  optimality. 

fiT.(  S 


Case  (b) : A . <h  , This  case  is  proved  in  a manner  similar  to  that 

aj  a— l,j 

used  for  case  (a)  above,  with  A . being  increased  from  0 to  x.  and 

a-l,j  J 

A , decreased  from  x.  to  0. 

»J  J 

Ey  repeated  use  of  cases  (a)  and  (b)  it  follows  that  a variable  A^ 
will  be  positive  in  an  optimal  solution  only  if  its  value  falls  between 

hi-l,j  *nd  hiJ  * 


f '*r*.*-. 


•■  vr  'SfV».’5i«r-.  '•'.rxT’ *%•—• 


' ' '•  «JU?>  '-|V*f  ; » 
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